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Abstract 

Theoretical models and rate equations relevant to the Soai reaction are reviewed. It is found that 
in a production of chiral molecules from an achiral substrate autocatalytic processes can induce 
either enantiomeric excess (ee) amplification or chiral symmetry breaking. Former terminology 
means that the final ee value is larger than the initial value but depends on this, whereas the 
latter means the selection of a unique value of the final ee, independent of the initial value. The 
ee amplification takes place in an irreversible reaction such that all the substrate molecules are 
converted to chiral products and the reaction comes to a halt. The chiral symmetry breaking is 
possible when recycling processes are incorporated. Reactions become reversible and the system 
relaxes slowly to a unique final state. The difference between the two behaviors is apparent in the 
flow diagram in the phase space of chiral molecule concentrations. The ee amplification takes place 
when the flow terminates on a line of fixed points (or a fixed line), whereas symmetry breaking 
corresponds to the dissolution of the fixed line accompanied by the appearance of fixed points. 
Relevance of the Soai reaction to the homochirality in life is also discussed. 

PACS numbers: 
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I. INTRODUCTION 



In the Human Genome Project (from the year 1990 to 2003) sequences of chemical base 
pairs that make up human DNA were intensively analyzed and determined, as it carries 
important genetic information. DNA is a polymer made up of a large number of deoxyri- 
bonucleotides, each of which composed of a nitrogenous base, a sugar and one or more 
phosphate groups Similar ribonucleotides polymerize to form RNA, which is also an 
important substance to produce a template for protein synthesis. RNA sometimes carries 
genetic information but rarely shows enzymatic functions [l]. 

One big issue of the post Genome Project is the proteomics since proteins play crucial 
roles in virtually all biological processes such as enzymatic catalysis, coordinated motion, 
mechanical support, etc. [lj]. A protein is a long polymer of amino acids, and folds into 
regular structures to show its biological function. 

Sugars and amino acids in these biological polymers contain carbon atoms, each of which 
is connected to four different groups, and consequently is able to take two kinds of stere- 
ostructures. The two stereoisomers are mirror images or enantiomers of each other, and 
are called the D- and L-isomers. Two isomers should have the same physical properties 
except optical responses to polarized light. Therefore, a simple symmetry argument leads 
to the conclusion that there are equal amount of D- and L-amino acids or sugars in life. But 
the standard biochemical textbook 1] tells us that "only L-amino acids are constituents of 

n 

oroteins" ( [lj], p. 17), and "nearly all naturally occurring sugars belong to the D-series" ( 
1|, p. 346). There is no explanation why the chiral symmetry is broken in life and how the 
homochirality has been brought about on earth. 

It is Pasteur who first recognized the chiral symmetry breaking in life in the middle of the 
19th century. By crystallizing optically inactive sodium anmonium racemates, he separated 
two enantiomers of sodium ammonium tartrates with opposite optical activities by means of 
their asymmetric crystalline shapes 2[ . Since the activity is observed even in the solution, it 
is concluded that the optical activity is due to the molecular asymmetry or chirality, not due 
to the crystalline symmetry. Because two enantiomers with different chiralities are identical 
in every chemical and physical properties except the optical activity, "artificial products 
have no molecular asymmetry" and in 1860 Pasteur stated that "the molecular asymmetry 
of natural organic products" establishes "the only well-marked line of demarcation that 
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can at present be drawn between the chemistry of dead matter and the chemistry of living 
matter" (3)]. Also by using the fact that asymmetric chemical agents react differently with 
two types of enantiomers, he separated enantiomers by fermentation. Therefore, once one 
has an asymmetric substance, further separation of two enantiomers or even a production 
of a single type of enantiomer may follow J4]. But how has the first asymmetric organic 
compound been chosen in a prebiotic world? 

This problem of the origin of homochirality in life has attracted attensions o 
scientists in relation to the origin of life itself since the discovery of Pasteur 
Japp expected a "directive force" when "life first arose" [sj], and various "directive forces" 
are proposed later such as different intensities of circularly polarized light in a primordial 
era, adsorption on optically active crystals, or the parity breaking in the weak interaction. 
However, the expected degree of chiral asymmetry or the value of the enantiomeric excess 
(ee) turns out to be very small 6|, and one needs a mechanism to amplify ee enormously to 
a level of homochirality. 

Another scenario for the origin of homochirality was suggested by Pearson 8| such that 
the chance breaks the chiral symmetry. Though the mean number of right- and left-handed 
enantiomers are the same, there is nonzero probability of deviation from the equal popu- 
lations of both enantiomers. The probability to establish homochirality in a macroscopic 



system is, of course, very small 



9|, but "chance produces a slight majority of one type of 



enantiomers and "asymmetric compounds when they have once arisen" act as "breeders, 
with a power of selecting of their own kind of asymmetry form" 8[. In this scenario, pro- 
duced enantiomer acts as a chiral catalyst for the production of its own kind and hence this 
process should be autocatalytic. 

Nearly a century later, a long sought autocatalytic system with spontaneous amplification 
of chiral asymmetry is found by Soai and his coworkers [h| : In a closed reactor with achiral 
substrates addition of a small amount of chiral products with a slight enantiomeric imbal- 
ance yields the final products with an overwhelmingly amplified ee [lO]. The autocatalytic 
system also shows significant ee amplification under a variety of organic and inorganic chiral 



initiators with a small enantiomeric imbalance 



ll| . When many reaction runs are performed 



without chiral additives, about half of the runs end up with the majority of one enantiomer, 
and the other half end up with the opposite enantiomer. The probability distribution of 
the ee is bimodal with double peaks, showing amplification of ee and further indicating the 
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occurrence of chiral symmetry breaking [l2, [jjj]. 

As for the theoretical research on the chiral symmetry breaking, Frank is the first to 
show that a linear autocatalysis with an antagonistic nonlinear chemical reaction can lead 
to homochirality [ijj]. His formulation with rate equations corresponds to the mean- field 
analysis of the phase transition in nonequilibrium situation 16[, and other variants have 



been proposed p. 



Q, Q, Q : 



20 
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22l | . All these analyses are carried out for open 
systems where the concentration of an achiral substrate is kept constant. Asymptotically 
;he system approaches a unique steady state (], 17, Q, 
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20| or even an oscilatory state 



22], independent of the initial condition. When the final state is chirally asymmetric, 



one may call this a chiral symmetry breaking in the sense of statistical mechanics, and it 
is indicated in the right upper corner in Table 1. On the other hand, the Soai reaction 
is performed in a closed system and an achiral substrate is converted to chiral products 
irreversibly: The substrate concentration decreases in time, and eventually the reaction 
comes to a halt. After the reaction the ee value increases but its final value is found to 
depend on the initial state. Even though the ee is amplified, the history dependent behavior 
is different from what we expect from phase transitions in statistical mechanics, so that 
we write down "ee amplification" in the left lower corner in Table 1. We have analyzed 
theoretically the chiral symmetry breaking in a closed system in general, and found that 
with an irreversible nonlinear autocatalysis the system may show ee amplification as in the 
Soai reaction. Furthermore, with an additional recycling back reaction, a unique final state 



with a finite ee value is found possible, which we ca 



in the right lower corner in Table 1 [2i 
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1 a chiral symmetry breaking as is shown 
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28lj. There are also many theoretical 
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Here we give a 



works on the Soai reaction in a closed system [29 
brief survey of various theoretical models of the ee amplification and the chiral symmetry 
breaking in a closed system, though it is by no means exhaustive. 



II. AMPLIFICATION IN MONOMERIC SYSTEMS 

Two stereostructural isomers are called D- and L-enantiomers for sugars and amino acids, 
but, for general organic compounds, R and S representation is in common. We adopt R — S 
representation hereafter. 

We consider a production of chiral enantiomers R and S from an achiral substrate A in a 
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TABLE I: The ee amplification and the chiral symmetry breaking in an open and a closed system. 





Initial Condition 


Dependent 


Independent 


System 


Open 




Chiral Symmetry Breaking 


Closed 


ee Amplification 


Chiral Symmetry Breaking 



closed system. Actually, in the Soai reaction, chiral molecules are produced by the reaction 
of two achiral reactants A and B a.s A + B — > R oy A + B — > S. But in a closed system a 
substrate of smaller amount controls the reaction, since it is first consumed up and holds the 
reaction to proceed further. Therefore, in order to grasp main features of chirality selection, 
it is sufficient to assume that an achiral substrate A of a smaller amount converts to R or S. 
The process may further involve formation of intermediate complexes, oligomers etc., but 
for the purpose to discuss about relations and difference between the amplification and the 
symmetry breaking in chirality, we restrict our consideration in this section to the simplest 
case where only monomers are involved. 



A. Reaction Schemes 



The spontaneous production of chiral molecules R or S from an achiral substrate A is 
described by reactions 

A H R, A^S. (1) 

We assume the same reaction rate ko for the nonautocatalytic spontaneous production of R 
and S enantiomers, since both of them are identical in every chemical and physical aspects: 
No advantage factor [6| is assumed. 

Since the spontaneous production ([1]) yields only a racemic mixture of two enantiomers, 
one has to assume some autocatalytic processes. The simplest is a linear autocatalysis with 
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a reaction coefficient k\ as 



A + R ^ 2R, A + S^2S. (2) 



There may be a linear cross catalytic reaction or erroneous linear catalysis with a coefficient 

as 

A + R^S + R, A + S^R + S. (3) 

As will be discussed later, the linear autocatalysis alone is insufficient not only to break the 
chiral symmetry but also to amplify ee. 

Beyond these linear autocatalyses, nonlinear effects such as quadratic autocatalysis have 



been considered as [6|, ]23j, 



and cross catalysis as 



30|; 



A + 2R^3R, A + 2S h 3S, (4) 



These hig 



A + R + S 2R + S, A + R + S ^ R + 2S. (5) 
rer order autocatalytic processes may actually be brought about by the dimer 



catalysts 181 ] . but we confine ourselves to the simplest description (jlj) and (]5]) in terms 
of only monomers R and S in this section. Consideration on dimers is postponed in the 
following sections. The reaction (j3J) alone can give rise to the amplification of chirality as 
will be discussed. One notes that all these processes ([!]) to (]5]) unidirectionally produce 
chiral enantiomers and thus the reaction comes to a halt when the whole achiral substrate 
is consumed; the total process is irreversible. 



B. Rate equations and an order parameter 

Concentrations r and s of two enantiomers R and S vary according to the reaction 
processes described in the previous subsection. In an open system, the achiral substrate is 
steadily supplied in such a way that its concentration a is kept constant. On the contrary, 
in a closed system, achiral substrates are transformed to chiral products and only the total 
concentration c = a + r + s of reactive chemical species is kept constant. By denoting the 
total concentration of monomeric chiral products with a subscript 1 as 

qi = r + s, (6) 
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the concentration of the substrate a is determined as 

a = c-q 1 , (7) 

which should be non-negative. Then, the development of reactions in a closed system can 
be depicted as a flow in the triangular region in the r-s phase space 

< r, s, r + s < c. (8) 

The reaction schemes (JTJ) to (J5J) induce rate equations for r and s as; 

dv 

— = [f{r) + k[s + k' 2 rs]a 

°± = [/(s ) + k > xr + k > 2rs]a (9 ) 

with an effective rate coefficient 

/(r) = k + hr + k 2 r 2 . (10) 

Since the right-hand side of Eq.Q for velocities f = dr/dt and s = ds/dt are always non- 
negative and in proportion to the concentration of the achiral substrate a, both r and s 
never decrease and stop to grow after a vanishes (a = 0). 

The asymptotic behavior of a two-dimensional autonomous dynamical system is, in gen- 
eral, known to have fixed points or lines where r = s = or to have limit cycles where 
f, s ^ 0. Since the system is irreversible as the concentrations of chiral products always 
increase at the cost of the substrate A as 

§ = -lf(r) + /(*) + K(r + s) + 2k' 2 rs]a < (11) 

and no epimerization process R <-> S is assumed, limit cycle is impossible: Production stops 
when all the substrate molecules are transformed to chiral ones; a(t = oo) = 0. Since 
r = s = for a = = c — r — s, the diagonal line r + s = c is a fixed line for this irreversible 
system. Any points on the fixed line are marginally stable because they have no guaranteed 
stability along the fixed line. 

Our interest is the dynamical behavior of ee or the chiral symmetry breaking and, in 
particular, as to whether ee increases asymptotically. In order to quantify the monomeric 
ee or the degree of chiral symmetry breaking, we define a monomeric order parameter as 

r — s r — s 

0i = — — = , 12 

r + s qi 



7 



which is zero for a symmetric or racemic state (r = s), and is nonzero for a symmetry-broken 
or chiral state. The evolution equation of <\>\ is immediately derived from Eq.Q as 



dt 



B^l (13) 



with coefficients 



A(t) = -2-(k + k[ qi ) + B(t), 
Qi 

B{t)= l -{k 2 -k' 2 ) qi a. (14) 

However, as a(t = oo) = 0, both the coefficients A and B vanish asymptotically, and one 
cannot determine the final value of ee from the evolution equation (fl3l) . On a fixed line 
r + s = c the ee value varies from 0i = — 1 at a point (r, s) = (0, c) to <pi = 1 at (c,0), and 
the final values of <fti can be anything between +1 and —1. In order to obtain the final value 
of ee, one has to solve the evolution (Q and especially to determine the trajectory of the 
dynamical flow in r — s phase space. 

C. Flow Trajectory in a r — s Phase Space 

As the first step to analyze flow trajectories in a phase space, we consider the sim- 
plest case when reactions such as spontaneous production (CQ), linearly autocatalytic (j2J) or 
quadratically autocatalytic (jlj) reactions are active, respectively. Then the rate equations 
are simplified as 

I = /(r) °' 

It = /(s)a ' (15) 

with the rate coefficient / defined in Eq.f llQI) . Flow trajectories are obtained by integrating 
a trajectory equation 

^ = M. (16) 
dr f(r) 1 ; 

In the following, an initial condition is set as r = vq, s = sq, ao = c — tq — sq (or q = q\$ = 
r + s , (f>i = 0i )O = (r — so)/<7i,0) a o = c ~ Qifi), an d we discuss how the final ee value 0i )OO 
depends on the initial value 01,0 for a few typical cases. 
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(a) 



(b) 



(c) 



FIG. 1: Flow diagrams for (a) spontaneous production (ko > 0, k\ = k 2 = 0), (b) linearly 
autocatalytic (fei > 0, &o = ^2 = 0), and (c) quadratically autocatalytic (&2 > 0, k$ = k\ = 0) 
reactions. 

1. Spontaneous Production 

When chiral molecules are produced only spontaneously and the system has no autocat- 
alytic processes (k > 0, k\ = k 2 = 0), the trajectory in the r — s phase space is obtained 
as 



The flow diagram is a straight line with a unit slope, terminating at the fixed line r + s = c, 
as shown in Figfl](a). Half of the substrate ao initially present convert to R and the 
remaining half to S, and the ee decreases to the value 0i jOO = (qi,o/ 'c)0i,o- The smaller the 
ratio qifi/c of the initial amount of chiral species to the total amount of active reactants is, 
the more the ee decreases, as shown in Figj2]^a). 

2. Linear Autocatalysis 

When only the linearly autocatalytic production of chiral molecules is active (fci > 0, k Q = 
k 2 = 0), the trajectory is obtained as 



r — s = 



r — s = const. 



(17) 




s s 



(18) 
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FIG. 2: The final ee (f>i >00 versus the initial ee 4>i,Q for (a) spontaneous production (ko > 0, ki = 
^2 = 0), (b) linearly autocatalytic (ki > 0, ko = ki = 0), and (c) quadratically autocatalytic 
(&2 > 0, ko = k± = 0) reactions. Each curve corresponds to a fixed ratio of the initial amount 
of chiral species to the total amount of active reactants qifi/c. As the ratio decreases, the initial 
slope of the curve decreases in (a), does not change in (b), and increases in (c). 

and the flow diagram is a straight line radiating from the origin and terminating at the fixed 
line r + s = c, as shown in Figfjjb). The ee value remains constant as 



'l.oo — Vi.o 



const. (19) 



and is independent of the initial ratio q^o/c, as shown in FigfS^b). 

If the erroneous linear catalysis (JSJ) is at work in addition to the linear autocatalysis (|2j), 
rate equations are modified as 

— = {hr + k x s)a, 
at 

^ = (jfe l5 + k[r)a. (20) 



These equations can be combined to yield 



dqi q± 

with a fidelity factor / defined by 



W)-, (21) 



= hJ< (22) 

J h + k[ y ' 
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which is unity without the error k[ = and vanishes when k\ = k[. The ee value 0i is easily 
integrated as 

01,oo = <£i,o (^7) < 01,0- (23) 

This result shows that the final value is less than the initial 0^0 whenever there is a finite 
error process / < 1. In particular, if / = 0, the final ee behaves the same with that of the 
spontaneous production. 



3. Quadratic Autocatalysis 

With only a quadratic autocatalysis (k% > 0, ko — k\ — 0), the trajectory is obtained as 

1111 

= = const. (24) 

r s r s 

and the flow diagram is hyperbolae radiating out from the origin and terminating at the 
fixed line r + s = c, as shown in FigfjTc). Below the diagonal r = s line, the flow bents 
downward and the ratio s/r decreases as time passes, and consequently the ee value 0i 
increases. In other words, by using the relation r = qi(l + 0i)/2 and s = q±(l — 0i)/2, the 
ee value is easily found from Eq.( |24l) to satisfy a relation 

Mi (25) 



1-0? 



with a constant A determined by an initial condition. For a given value of qi, the ee value 0i 
is determined graphically as a cross point of two curves representing both sides of Eq. (l25l) . 
as shown in Figj3j Then, it is evident that <pi increases as qi increases. The asymptotic 
value corresponding to q± = c is obtained as 

1 



>l,oo 



2cA 



1 + (2cA) 2 - 1 (26) 



The relation between the initial value 0i o of the ee and its final value 0i jOO is depicted 
in Figj2](c) for fixed values of the initial amount g^o of chiral species relative to the total 
amount c. The smaller the initial ratio qi,o/c, the more prominent the ee amplification. 

The ee amplification curves possess a remarkable property such that they depend on the 
initial ratio qi^/c but are independent of the rate coefficient k 2 . The latter determines the 
time evolution of <pi but not its final value. If there are other processes involved, for instance 
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FIG. 3: The ee value (j>i is given at a cross point with two curves, y = 4>i/(l — 0f) and y = Aq\. 

the spontaneous production ko 7^ 0, then amplification curves should depend on the ratio 
k /k 2 c 2 . 

If there is a racemizing cross catalysis (jSJ) with k' 2 , the trajectory equation is rewritten as 



dqi fx 1 + <i>l 91 

with a crossing effect factor defined by 



(27) 



= l^t <28) 

which is unity without a cross catalysis term [k! 2 = 0) and vanishes when k 2 = ^2- The 
solution of Eq. f[2"T|) is given by 

(1 _ (fiyi+f* )/2 = Ag i x • ( 29 ) 

From this result one can see that the asymptotic value 0i jOO decreases as f x decreases. In 
particular, when f x — 0, is independent of gi and does not vary; 0i jOO = 0i )O . Namely, 
chirality amplification by k 2 is completely cancelled out by the racemization effect caused 
by k' 2 . 

III. AMPLIFICATION WITH HOMODIMER CATALYST 

In order to understand molecular mechanisms how the quadratic autocatalysis is brought 
about, the concept of dimer catalyst introduced by Kagan and coworkers may be relevant 
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32] . Assume that monomers R and S react to form homodimers R 2 and S 2 with 



a rate u\ loni and a heterodimer RS with a rate i^ et 

R + R"^R 2 , S + S u ^S 2 , R + S^RS. (30) 
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Corresponding decomposition processes are described as 

R 2 "h>R + R, s 2 ^S + S, RS^R + S. (31) 

Homodimers thus formed are assumed to catalyze of their own enantiomeric type as 

A + R 2 fc ^ R + R 2 , A + S 2 k ±T S + S 2 , (32) 

whereas heterodimers have no preference in enantioselectivity as 

A + RS ^ R + RS, A + RS ^ S + RS (33) 

Now, the state of the system is described not only by the concentrations of monomers 
r and s, but also by those of dimers denoted as [R 2 ], [S 2 ] and [RS] respectively. The rate 
equations are 

dv 

— = k a + k hom a[R 2 ] - 2u hom r 2 + 2v' hom [R 2 ] 



dt 
d[RS] 



vnetrs + is' hct [RS] + k het a[RS], 
Vhomr 2 — ^[^s], 

u het rs - i/ het [RS], (34) 



dt 

and the corresponding ones for the S enantiomer. The ee or the chiral order parameter is 
now defined by 

, r + 2[R 2 ] - s - 2[S 2 ] 



q 

where q is the total concentration of chiral molecules given as 



(35) 



q = r + s + 2[R 2 ]+2[RS}+2[S 2 ). (36) 

If the dimerization and decomposition proceed very fast, the dimer concentrations take 
quasi-equilibrium values which satisfy [R 2 ] = [RS] = [S 2 ] = and are determined by the 
instantaneous monomer concentrations r and s as 

r n 1 ^hom 2 rnpi ^het r a i ^hom 2 fon\ 

[R 2 ] = — — r , [RS] = — rs, [S 2 ] = — — s . (37) 

^hom ^het ^hom 

Then, the equations for monomers in Eq.( l34|) reduce to Eq.flH]) with the coefficients k 2 and 
k' 2 given as 

k 2 = ^hom - ; — j k 2 = khct— — (38) 
1/1/ 

horn het 
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and with k\ = k[ = 0. 

As the reaction proceeds, the whole achiral substrate is ultimately transformed to chiral 
products so that a = c — = 0. This represents a quadratic curve of fixed points in r — s 
phase space. When there is solely a homodimer catalytic effect k 2 > without spontaneous 
production nor heterodimer racemization k = k' 2 = 0, the initial state with (r ,So) follows 
the hyperbolic flow trajectory (124)) and ends up on the curve of fixed points, q = c, as shown 
in Fig Ufa). 

The ratio r/s increases along the flow trajectory if it is initially larger than unity. Thus 
the r — s asymmetry increases, and the final value of the ee \4>oo\ is larger than the initial 
value 1 0o I- Figure H£b) shows 0o ~~ <Poo relation for cases of positive <fi. The enhancement is 
more evident when the initial ratio q^/c is smaller. These essential features are the same 
with the monomeric system. 




FIG. 4: (a) Flow diagrams for homodimer catalyst with iWnC/^hom = 0-1 an d ^et = ^hct = 0- 
Flow terminates on a curve of fixed points q = c. (b) The final ee value <j)oo as a function of its 
initial value <po for various amount of initial chiral species qo. The less qo is, the more amplified 
0oo is. 



IV. AMPLIFICATION WITH ANTAGONISTIC HETERODIMER 



Another model for the chiral asymmetry amplification is a variant of the Frank model 
in a closed system 33(. In its simplest form, the model consists of reaction schemes of a 
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spontaneous production (JTJ, a linear autocatalysis ([2]) and a heterodimer formation (I3T1) . 
Then the rate equations for monomers are the same with those of the original Frank model 
[is] ], but they are supplemented with the heterodimer formation 

dr 



— - = k a + k x ar - u he trs, 
dt 



ds 

~dt 
d[RS] 

dt 



k a + kias — z/ hct rs 



(39) 




,'r = s 



0.2 0.4 0.6 0.8 1 

0o 



(a) 



(b) 



FIG. 5: (a) 3D flow diagrams for reactions with antagonistic heterodimers RS. Red lines represent 
fixed lines. Initial amount of achiral substrate is 0.2c, and the remaining qo = 0.8c is shared by 
R and S monomers. There are no heterodimers initially. Rate coefficients are ko = 0, k\c = 1, 
u hct c = 10. (b) Final 66 4>oo versus its initial value (fro for a given initial ratio qq / c. As Qo/c 
decreases, the ee amplification increases. 

Since every achiral substrate is consumed up eventually a(t = 00) = and all the reactions 
stop asymptotically, Eq. (1391) tells us that the product rs should vanish. If there is more R 
than S initially, S monomer disappears ultimately, for instance. But 5* molecules do not 
disappear nor decomposed back into achiral substrate. They are only incorporated into the 
heterodimer RS. The system is not determined solely by monomer concentrations r and s, 
(or <fri and q\) but it depends on heterodimer concentration [RS] as well. The flow takes 
place in a three-dimensional phase space of r, s, [RS], as shown in Fig[5ta). 
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The asymptotics discussed above forces both the concentrations of the substrate a = 
c — r — s — 2[RS] and the product rs vanish ultimately. These two conditions define a 
line of fixed points in the three-dimensional r — s — [RS] phase space. If the initial state 
has a prejudice to R enantiomer such as r > s , then the system ends up on a fixed line 
r + 2[-RjS] = c on a s = plane, as shown in FigJS^a). Otherwise with r < Sq, the system 
flows to another fixed line s + 2[i?S"] = c on a r = plane. The plane r = s is a dividing 
surface that determines the final preference of chirality. 

One can readily find that ee is amplified by the reaction involving antagonistic het- 
erodimer, as shown in Figj5](b). For a given amount of initial total concentration of chiral 
species q = r + s , the ee value increases by the completion of a reaction. No het- 
erodimer is assumed to exist initially, [RS]q = 0. The amplification is more significant for 
smaller amount of qo/c and 0o- We note that the rate constant of heterodimerization 
has to be large in order to produce strong ee amplification. Minority enantiomer should be 
quickly incorporated into heterodimers before the achiral substrate is consumed up. Other- 
wise, both enantiomers consume up the achiral resource A and slowly forms heterodimers 
afterwards but without significant ee amplification. For a large z/h et , minority enentiomers 
are incorporated into heterodimers in an early stage. Therefore, flow trajectories first point 
upwards in Figj5]^a), indicating the rapid increase of heterodimer concentration [RS]. The 
trajectory is in fact inclined such that the minority enantiomer decreases its monomeric 
concentration. When it almost disappears, then the trajectory makes a quick turn in the 
vicinity of r = plane or s = plane, indicating the concentration growth of the surviving 
majority enantiomer. 

V. RECYCLING AND FLOW: CHIRAL SYMMETRY BREAKING 

So far, we have discussed several models for the ee amplification. There, the irreversibility 
of the dynamical system is essential to make the final state dependent on the initial condition. 
In this section, we consider the effect induced by introducing reversibility in the system by 
including various types of recycling processes. For simplicity, we consider in terms of a 
monomeric system. 

In addition to the catalytic production processes ([T]) to (j5J), we include the recycling 
process such that enantiomeric products R or S undergo a back reaction to the substrate A 
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(a) (b) (c) 

FIG. 6: 4>\ versus the velocity 0i = A4>\ — B(f>\ for three typical cases: (a) A < B < 0, (b) 
A < < B , and (c) < A < B. 



231 ] with a reaction rate A as 



R^A, S ^ A. (40) 

We call this reaction a linear recycling hereafter. One can also assume a nonlinear type 
of recycling process such that different chiral species by their encounter react back to 
achiral molecules as 

R + S^2A. (41) 
This back reaction can be a combined effect of the following two processes 

R + S ^ R + A 

R + S A S + A. (42) 

The rate equations now read as 

= [/( r ) + K s + k' 2 rs]a — firs — \r 

ds 

— = [f(s) + k[r + k' 2 rs]a — firs — As. (43) 

(JjO 



The term with fi has a same form with the contribution of Frank's mutual inhibition [15| . The 
corresponding time-evolution of <pi is readily derived in the form of Eq. f|T3|) with coefficients 
A and B determined by symmetric quantities q± and a = c — q\ as 

A(t) = -2-(k + k[ qi ) + B(t), 

B(t) = ^[(k 2 -k' 2 )a + fx] qi . (44) 
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The evolution E g.(fT3l) of the order parameter has a similar form with the time-dependent 
Landau equation 16j which is fundamental in nonequilibrium phase transitions. The asymp- 
totic value of the order parameter 0i iOO is determined as the zero of the velocity 4>i- The 
main difference from the standard model of phase transitions lies in the time- dependence 
in the coefficients A(t) and B(t) induced by that of the achiral concentration a(t) and the 
total chiral concentration qi(t). Because the concentrations a and q± are non-negative, A(t) 
cannot exceed B(t); A(t) < B(t). 

In previous sections we have focused our studies on the cases with A = [/, = 0. In these 
cases, the asymptotic value of a vanishes, a(t = oo) = so that A(t = oo) = B(t = oo) = 0, 
and the asymptotic value of the order parameter 0i cannot be determined from Eq. (jl3j) . 
Another example with A = B = happens for a special case with ko = k[ = fi = and 
/c2 = k' 2 . In this case, we can calculate trajectories and find a fixed line shifted from a = 
to the one with a finite value of a(t = oo) 27j|: Flow trajectories are along lines passing 
through the origin, similar to those shown in Fig{T^b), but terminating at a shifted fixed 
line. Since a = is no more a fixed line, trajectories not only point out from the origin, but 
also come out from the line a = 0. 

In all the other cases, with either a linear or a nonlinear recycling process or with both, 
the coefficients A and B are no longer zero at the same time, and a definite value of the 
order parameter 0i is obtained asymptotically. The reason for the above is as follows. If 
the nonlinear recycling exists as \x > 0, B becomes nonzero since fJ,qi }0 o > 0. If the linear 
recycling exists as A > 0, not all the achiral substrate transform to chiral products but a 
finite amount remains asymptotically as a(t = oo) > 0. Therefore, nonzero values of ko, k[ 
or k 2 7^ k' 2 give contributions to the coefficients A or B. 

For nonzero A and/or B, the velocity <pi = d<pi/dt has three typical behaviors as a function 
of the order parameter 1; as shown in Figj6j In the region with a positive (negative) velocity 
0i, 0i increases (decreases) and moves to the right (left), as indicated by arrows on the 0i 
axis in FigEl When A is negative, as in Figj6]^a) or (b), 0i eventually approaches the racemic 
value 0i = within its range of definition |0i| < 1: The racemic state is a stable fixed point 



for A < 0. Equation (jUj) tells us that this happens when the spontaneous production fc 
or error in a linear autocatalytic process k[ are sufficiently strong, or quadratic catalysis k 2 
is weaker than the cross catalysis k' 2 . The coefficient of the linear autocatalysis k\ as well 
as the linear recycling process A are absent in the coefficients A and B, and cannot affect 
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directly the chirality of the system. These linear processes with k\ and A affect the chirality 
only implicitly through a and q\. 

When A is positive, as in the case of FigJ^c), the coefficient of the cubic term B is also 
positive, and the velocity <f) vanishes at three values of <pi in the range of —1 < (pi < 1. 
This is possible if a strong quadratic autocatalysis &2 > k' 2 exists together with a linear 
recycling A > 0, or if a linear autocatalysis k\ and a nonlinear recycling /i > coexist. By 
following the direction indicated by the arrows for positive <pi, the order parameter ends up 
at a definite value 



Mt = oo) 

n ^ = "W^) (45) 

and the chiral symmetry breaking takes place. If the system starts with a negative 0i, it stops 
at — 0i,oo- The final state of the system depends on the initial sign of 01,0, but independent 
of its magnitude |0i,o|- We note that the asymptotic values of the coefficients A(t = oo) and 
B(t = oo) are independent of the initial condition and so is |0i )O o|- This situation is what 
we mean the "chiral symmetry breaking" in the sense of statistical mechanics. 

We first discuss the chiral symmetry breaking when there is a nonlinear recycling, /i > 
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27| . The rate Eq. (T4l?j) tells us that a = c — r — s = is no more a fixed line, but we 
have many fixed points. The simplest is the case with only a linear autocatalysis k\ > but 
nothing else (k = k[ = k 2 = k' 2 = A = 0). In this case we have four fixed points in the r — s 
phase space: two unstable ones, O and U, and two stable ones at (r, s) = (c, 0) or (0, c), as 
is shown in FigJT^a). At the stable fixed points the final ee take values 0i iO o = ±1- These 
values agree with those expected from Eq. p5j) with A(t) = B(t) = fiqi/2 as determined 
from Eq. fl44p . The nonlinear recycling (42) returns the same amount of enantiomers back to 
the achiral substrate, but the relative decrement in concentration is larger for the minority 
enantiomer than for the majority one. The damage is further amplified by a large reduction 
of linear autocatalytic effect for the minority enantiomer, and it disappears at last. 

As for the second example, we consider the case with a quadratic autocatalysis /C2 > and 
a linear recycling A > 0. Because the linear recycling takes place whenever there is nonzero 
concentrations of enantiomers, a few achiral substrate always remain, a(t = oo) ^ 0. Thus 
the diagonal line q\ = r + s = c is no more a fixed line. Instead there appear fixed points 
in general. In the present case with a finite A > 0, totally seven fixed points appear; three 
stable (O, Si 2 ) and four unstable ones (Ui_ 4 ), as shown in FigJTl As reaction proceeds, the 
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FIG. 7: Flow diagrams with (a) a linear autocatalysis and a nonlinear recycling ( k\ = [i > 0, ko = 
&2 = A = ) similar to the Frank model, and (b) a quadratic autocatalysis and a linear recycling 
(k 2 c 2 = A > 0, k = kx = n = 0). 

system approaches a state associated with one of the stable fixed points. The origin O is 
not interesting, since all the chiral products are recycled back to the achiral substrate. Also 
its influence extends only in a small region around O, bounded by Ui_3, if the life time 1/A 
of chiral products is long enough. Close to this fixed point, influence of the spontaneous 
production process with k should be taken into account, anyhow. 

When the initial system is sufficiently away from O, the system approaches Si or S2; Both 
correspond to homochiral states with 0i jOO = ±1. When the initial configuration is close to 
the racemic state or the diagonal line r — s, the system approaches a racemic fixed point 
U4 at first. But, while the recycling process returns the chiral enantiomers back to achiral 
substrate, the majority enantiomer increases its population at the cost of the minority one 
along the flow curve diverging from U 4 to Si or to S 2 , and eventually the whole system 
becomes homochiral. 

One may wonder whether recycling processes such as a linear or a nonlinear back reaction 
exist in relevant autocatalytic systems. So far, we are not aware of their existence. It is, 
however, possible that the back reaction rates A or fi are nonzero but exceedingly small to 
be detected in laboratory experiments. Concerning with the problem of homochirality in 
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life, very small A or /i are not unimaginable, considering the geological time scale for its 
establishment on earth. 

On the other hand, a possibility to provide a system with a recycling process is proposed 



theoretically 26j: One simply let the reaction with ee amplification run in an open flow 
reactor. In an open system there is reactant flows between the system and the environment 
E. Let the solution with the achiral substrate A be supplied by an inflow with an influx F 

as 

E ^ A (46) 

and the products R and S as well as the substrate A be taken out by an outflow with a flow 
out rate A as 

A^E, R^E, S^E. (47) 
In terms of rate equations the processes are described as 

= f(r)a — Ar, ^£ _ f( s \ a _ \ s 
dt y J dt v ' 

^ = -[f(r) + f(s)]a + F-\a. (48) 

Since the total concentration a + r + s follows the time evolution d(a + r + s)/dt = F — 
A(a + r + s), it approaches the steady state value F/X with a relaxation time 1/A. This is 
a consequence of unbiased outflow ( j48|) of all reactants with the same rate A. Consequently, 
even though we are dealing with an open system under a flow, the analysis is similar to 
the closed system by replacing the total concentration c with the steady state value F/X. 
Instead of recycling, therefore, constant supply of the substrate allows the system to reach a 
certain fixed point with a definite value of the order parameter 1; independent of the initial 
condition. 



VI. SUMMARY AND DISCUSSIONS 

We have presented various simple scenarios that we are aware of in relevance to the 
ee amplification of the Soai reaction; a quadratic autocatalytic model in a monomer or a 
homodimer system, or a linear autocatalytic model in an antagonistic heterodimer system. 
All these models can realize ee amplification such that the final value of the ee \4>cx>\ depends 
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on but is larger than the initial value \4>o\, as schematically shown in FigJBJ The curve in 
the figure represents 0o — 0oo relation for a given initial ratio qo/c, namely the ratio of the 
amount go of the total chiral initiators R and S relative to that of the total reactants c. 
Amplification is more enhanced if the ratio of the chiral initiator q /c is smaller. This plot 
also shows the possibility to increase the final ee by repeating the reaction. 

To achieve high ee values, Soai et al. performed experiments of consecutive autocatalytic 
reactions, where after each run of the reactions, resulting solutions are quenched by adding 
acid and purified and reused in the next reactions. How the high ee value is attained in 
this type of experimental procedures can be explained illustratively in Fig. 8. By choosing 
a specified initial condition 4>^\ the autocatalytic reaction ends up giving a final ee as 
4> — 4>oo determined by the intersection between the vertical line starting from point 
and the round curve. The consecutive reaction starts with the purified products with the ee 
as the initial ee . This feature is illustrated by extending a horizontal line from <j)2o on 

(2) 

the round curve to the 0o = 4>oo straight slope. This second reaction yields the ee value 0oo , 
which is determined by the intersection between the vertical line starting from (p^ = 
point and the round curve. To repeat these procedures, the ee approaches rapidly the perfect 
|0oc°' ) | = 1 with the speed depending on the curvature of the curve, which is determined by 
qo/c. This is a mechanism of ee amplification and consecutive autocatalytic reactions, in 
contrast with the chiral symmetry breaking in the sense of statistical mechanics, where the 
final ee is unique and is reached by a single procedure. 

In all the cases considered, autocatalytic processes must be present, whether it be linear 
or nonlinear. To understand the actual mechanism of autocatalysis for the Soai reaction, 
the identification of the process at a molecular level is necessary, but that is out of scope of 
the present review. 

We would like to point out another interesting feature, namely the spatial aspect of the 



reaction 
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28 



36J. In usual experiments, the reactor is stirred to keep the system ho- 



mogeneous. If there is no stirring, inhomogeneous patterns may develop between domains 
of two enantiomers. Analysis of these patterns may give some insight into the mechanism 
of reaction dynamics. Let us imagine, for instance, a time evolution of patterns of chiral 
domains for diffusionless systems with ee amplification and with a chiral symmetry breaking. 
In an initial stage chiral domains are nucleated randomly, and the domain pattern may look 
similar for both systems. After domains come in touch with each other, they behave dif- 
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FIG. 8: Schematics of enantiomeric excess amplification. 

ferently: For a system with ee amplification, domain pattern freezes because there remains 
no achiral substrate for further reactions. On the other hand, for a system with the chiral 
symmetry breaking, neighboring domains compete with each other at domain boundaries to 
extend their own territories, and the domain pattern evolves to selection slowly. When dif- 
fusion or hydrodynamic flow effects are taken into account, the system with ee amplification 
simply homogenizes and the ee value shall be averaged out. For the system with the chiral 
symmetry breaking, these transport processes speed up the chirality selection and enhance 



the final value of ee 



28 



36]. 



So far we have discussed the ee amplification by starting from an initial state with a finite 
enantiomeric imbalance, and the effect of spontaneous reaction is mostly neglected, ko = 0. 
Towever, reaction experiments without adding chiral substances do produce chiral species 
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141 ] . implying that the coefficient ko of the spontaneous reaction should be nonzero. 
In these experiments, it is found that values of the ee order parameter of many runs are 
distributed widely between —1 and +1, and the probability distribution function has double 
peaks at positive and negative values of 0. 

There is a theoretical study on the asymptotic shape of probability distribution for nonau- 
tocatalytic and linearly autocatalytic systems with a specific initial condition of no chiral 



enantiomers 
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35| . Even though no ee amplification is expected in these cases, the proba- 
bility distribution with a linear autocatalysis has symmetric double peaks at <fi = ±1 when 
ko is far smaller than ki, ko <C k\/N with N being the total number of all reactive chem- 
ical species, A, R and S. This can be explained by the "single-mother" scenario for the 
realization of homo chirality, as follows. From a completely achiral state, one of the chiral 
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molecule, say R, is produced spontaneously and randomly after an average time l/2koN. 
Then, the second R is produced by the autocatalytic process whereas for the production 
of the first S molecule the spontaneous production is necessary. The waiting time for the 
second R is shorter than that for the first S by a factor about k^/ki . If this factor is far 
smaller than 1/N, there is only a negligible chance for the production of S enantiomer until 
all the achiral substrates turn into R. All the chiral products are the descendants of the first 
mother R produced by the autocatalysis, and no second mother is born by the spontaneous 
production: This is a single-mother scenario. It strongly reflects the stochastic character of 
the chemical reaction, and is not contained in the average description using rate equations. 
Further studies are necessary as for the stochastic aspects of the Soai reaction 37]. 

In relation to the origin of homochirality in life, an interesting question is whether the 
ee amplification is sufficient to achieve homochirality. From the experimental results of the 
Soai reaction without chiral initiator, the magnitude of the final ee \<poo\ is distributed from 
a small value to a value close to unity. It is rather difficult to imagine that a value close to 
unity has been accidentally selected on various places on earth in the prebiotic era. It seems 
rather natural that chiral symmetry breaking took place and a unique value of |0| has been 
finally selected. As has been discussed in the previous section, the flow in the open system 
can alter the system with ee amplification to the one with the chiral symmetry breaking. 
So we may say that the system with ee amplification in an open flow is one possibility. 



Recently, there appears several theoretical proposals such as polymerization models 26, 



38, 



391 ] or a polymerization-epimerization model 40[ to realize the chirality selection in polymer 



synthesis. Although we are still far away from satisfactory understanding of the origin 
of homochirality in life, the discovery of the autocatalytic reaction by Soai gives us great 
impetus and a sense of reality to this problem. 
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